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Abstract 

Recent claims in the literature have suggested that the WMAP quadrupole is not primordial 
in origin, and arises from an aliasing of the much larger dipole field because of incorrect satellite 
pointing. We attempt to reproduce this result and delineate the key physics leading to the effect. 
We find that, even if real, the induced quadrupole would be smaller than claimed. We discuss 
reasons why the WMAP data are unlikely to suffer from this particular systematic effect, including 
the implications for observations of point sources. Given this evidence against the reality of the 
effect, the similarity between the pointing-offset-induced signal and the actual quadrupole then 
appears to be quite puzzling. However, we find that the effect arises from a convolution between 
the gradient of the dipole field and anisotropic coverage of the scan direction at each pixel. There 
is something of a directional conspiracy here - the dipole signal lies close to the Ecliptic Plane, 
and its direction, together with the WMAP scan strategy, results in a strong coupling to the l2,-i 
component in Ecliptic co-ordinates. The dominant strength of this component in the measured 
quadrupole suggests that one should exercise increased caution in interpreting its estimated am- 
plitude. The Planck satellite has a different scan strategy which does not so directly couple the 
dipole and quadrupole in this way and will soon provide an independent measurement. 
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I. INTRODUCTION 



Measurements of Cosmic Microwave Background (CMB) anisotropics from space have 
laid the foundations of the standard model of cosmology. These observations provide prima 
facie evidence that the Universe is close to spatially flat with nearly scale-invariant initial 
density fluctuations. It is remarkable that only a handful of other parameters, specifying 
the fraction of baryonic and dark matter components, together with the local expansion rate 
and the redshift of reionization, can flt the data over such a wide range of scales. 

Despite the success of the standard model (see e.g cosmologists have searched for 
evidence of new physics. Unfortunately this practice is fraught with uncertainty because of 
the prevalence of a posteriori statistics. One particular avenue of investigation has focused 
on the CMB anisotropics at large angular scales. A low quadrupole, for example, was 
flrst noted by the Cosmic Background Explorer (COBE) m and subsequently conflrmed 
by the Wilkinson Microwave Anisotropy Probe (WMAP) j3|. Several other features have 
been discovered in WMAP data (see |j] and references therein), with subsequent debate 
about their statistical signiflcance. Hence it is crucial to investigate all potential sources of 
systematic error which could be important at large angular scales. 

Recently, it has been claimed that an important eflFect has been overlooked in the WMAP 
analysis |5|]. This particular issue relates to a 25.6 ms offset between recording the pointing 
and differential temperature data from the satellite, which translates into an angular error 
of about 7^ As part of the WMAP processing pipeline a dipole signal is flrst removed from 
the time-ordered data (TOD), and this has a much higher amplitude than the primordial 
anisotropics. If the pointing is incorrect a residual signal will remain in the TOD, which has 
been found to induce a quadrupole pattern in the flnal temperature sky maps. Interestingly, 
this quadrupole has similar spherical harmonic coefficients to those of the primordial 
signal measured by WMAP. 

In this short article we show why the claimed result cannot be correct. In doing so 
we uncover the physics behind the coupling of such systematic effects to certain harmonic 
modes. We also investigate the implications for the Planck experiment, which measures 
absolute rather than differential data, and scans the sky very differently to WMAP. 

II. WMAP RESULTS 

The WMAP TOD vector d can be written as 



where M is the mapping matrix, t the sky-map and n the radiometer noise [y]. Ignoring 
the noise, the outputs of the two radiometers are [3] 



where i is the Stokes / parameter and the pointing on the sky is denoted by fi, with the 
subscripts indicating the A and B-side beams. The pre-factor (1 ±Xim) for each term is due 
to transmission imbalance between radiometers, the latest values for which can be found 
in jsl. These factors are constant in time, and are easily absorbed into the mapping matrix. 

For the purposes of this work we assume the noise is zero (including noise would only 
effect the statistics of the results, not the general features). The maximum likelihood sky- 
map solution is then t = S to, where S = (M^ M) ^ and the 'iteration-0' map to = d. 



d = Mt + n 



(1) 



d = (1 + Xim) i(nA) - (1 - ^im) i(nB) 



(2) 
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FIG. 1: (Top): Induced quadrupole from a 25.6ms pointing offset in Galactic (left) and Ecliptic 
(right) co-ordinates. (Bottom): WMAP quadrupole derived from the 7- year ILC map in Galactic 
(left) and Ecliptic (right) co-ordinates. Units are in mK for each map. The bottom panels are 
strongly correlated with the upper ones, although with an amplitude about 3 times higher. 

The mapping matrix M has A^t rows and A^p columns, where A^t is the number of time 
samples for each radiometer and A^p the number of pixels in the map. 

To study the effects of a pointing offset, we assume our TOD contains a pure dipole 
signal. We fix the barycentric component to an amplitude 3.3463 mK with Galactic co- 
ordinates (/,6) = (263.87°, 48.2°) (g). With a CMB monopole temperature of Tq = 2.725 K, 
this corresponds to a Galactic velocity Vgai = (—26.2,244.0,274.4) kms~^. We also include 
an additional varying dipole component due the satellite motion around the Sun, such that 
the total velocity is v = Vgai + Vgoi. The dipole at each pointing is then (Tq/c) v • fi in 
temperature units. 

From this simulated time-stream we then remove the dipole evaluated at a small pointing 
offset (corresponding to the claimed timing discrepancy), leaving a residual signal in the 
TOD. We next run our map-making code, solving for the sky-map by a conjugate gradient 
method. We produce maps at a HEALPIX [lo| resolution of nside = 16. 

Maps at W-band, which contain 30 observations (each lasting 51.2 ms) per 1.536 second 
'science frame', are shown in Fig. [H In this case the central point of each observation 
corresponds to the claimed offset between recording the pointing information and recording 
the differential data from the satellite [5]. For W-band we can therefore use available WMAP 
software [l5| to compute the pointing offset, as the code also gives the interpolated pointing 
at the centre of each observation. The resulting map is characterized by a quadrupole with an 
amplitude ^ 6.5 /xK, and in Ecliptic coordinates is close to a pure I2, -i spherical harmonic 
mode. For comparison, we also plot in Fig. [1] (bottom panels) the WMAP quadrupole 
derived from the 7-year internal linear combination (ILC) map. 
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FIG. 2: Modulus of the scan tangent vector expectation value for WMAP (left), along with the 
underlying dipole field (right, units are mK), in Ecliptic co-ordinates. 

The induced quadrupole maps are similar in structure to the published WMAP 
quadrupole. However, we note that we obtain an amplitude lower than that claimed in j5|, 
by approximately a factor of 3. Hence, even if the claimed effect was real it could not have 
a large enough amplitude to explain the WMAP result. Nevertheless, the fact that it would 
produce a quadrupole so similar is intriguing and seems worthy of further investigation. 



III. UNDERSTANDING THE INDUCED QUADRUPOLE 

To make sense of this effect we built a toy model of the scanning strategy for satellites 
such as WMAP and Planck. The pointing vector can be obtained as a function of time by 
orienting the satellite focal plane (FP) in the +Zfp direction. The beam (or pair of beams in 
the case of WMAP) is specified by the opening angle with respect to the +^fp direction. 
In this frame the satellite spins at an angular velocity uj^ = 2it /t^^ where tg is the spin period. 

In order to achieve coverage at the Ecliptic poles the focal plane is rotated at a precession 
angle with respect to the anti-solar direction (which is now defined as the new spacecraft 
(SC) +^sc direction). The satellite then spins at a rate = 2ti /t^ in this frame. Finally, 
these co-ordinates are rotated into the Ecliptic Plane (with -\-z corresponding to the north 
Ecliptic pole) and the solar rotation rate is fixed at w^. 

This toy model accurately reproduces the pointing of WMAP and Planck. For WMAP 
311, the relevant parameters are = 70°, 6^ = 22.5°, tg = 129 s, tp = 3600 s, while for Planck 
ll\ they are ^b = 85°, 9p = 7.5°, = 60s and tp = 6 months. Note that the precession 
period is much longer for Planck than WMAP. 

We can introduce a time lag St to characterize the pointing offset - the residual signal 
at each pointing is then (Tq/c) 5tv • dn/dt. For a single beam experiment such as Planck 
the zeroth order approximation for the maximum likelihood solution map is the expectation 
value (in the long-time limit) of the scan tangent vector at each pixel (which we define as 
(dn/dt)) dotted with the dipole velocity. The overall amplitude is proportional to the spin 
rate of the satellite. 

For an experiment such as WMAP the to map is constructed from a difference signal 
between the two beams, {To/c)6tv • [driA/dt — driB/dt] (ignoring transmission imbalance 
factors). The maximum likelihood sky-map solution t is also complicated by the fact that 
the S matrix is not diagonal. We tested the implications of this by assuming that WMAP 
could also operate as a total intensity experiment. 



4 



In the default 'differencing mode' with the WMAP scan pattern and St = 25.6 ms, we 
recover the induced quadrupole shown in Fig. [1] to good precision. In 'total intensity mode' 
the induced signal was around a factor of several lower. In this case there is a subdominant 
contribution from the _i harmonic, which has an amplitude of around a factor of 4 less 
than the I2, -1 mode. The reason for these results is that, for WMAP^ dxix/dt o:^ —dii^/dt^ 
so the to map is similar in both differencing and total intensity modes. The to map contains 

-1 and I4, -1 components, but the latter is suppressed by the map-making operation in 
the differencing mode. 

The overall structure of the final t maps can therefore be understood by considering 
{dn/dt) and the underlying barycentric dipole field. We show these quantities in Ecliptic 
co-ordinates in Fig. [21 which should be compared with the induced map in Fig. [H Near 
the Ecliptic poles pixel scans are approximately isotropic, so the sum of dot products with 
the dipole velocity cancel over the course of a year's observations. Along the Plane, the 
component of the dipole velocity parallel to the scan direction is small, since pixels are 
mainly scanned orthogonal to the plane. Moreover, the sum of north/south scans cancel 
over a year. For angles at around 9^ from the poles pixels are mainly scanned in a single 
direction parallel to the equator. Here the dot product with the dipole velocity is large in 
certain azimuthal directions, resulting in a large induced signal. 

We note that it is an unfortunate coincidence that the direction of the dipole field in the 
Ecliptic Plane results in a quadrupole with approximately the same phase as the primordial 
signal of WMAP. For the WMAP scan strategy the dipole and quadrupole patterns therefore 
have a special and unfortunate relationship. Although the pointing offset effect appears not 
to be real, one can think of other systematics which could potentially couple to the large 
dipole field in a similar way. These include an asymmetric beam or a time constant in the 
instrument response. Therefore, it seems reasonable that one should take additional care 
in removal of the dipole from the TOD, since the induced signal would be similar to the 
primordial quadrupole. 

With our toy model in hand, the scan parameters can be adjusted to test the impact 
of other scan strategies on the induced quadrupole map. For the approach of the Planck 
satellite the pixel hits map is no longer azimuthally symmetric (after a year, say) due to 
the long precession period. The modulus of the scan tangent vector, along with the induced 
signal for a time delay of 6t = 25.6 ms, are shown in Fig. [31 The amplitude of the map is 
similar to that obtained for WMAP^ but the structure is very different due to the different 
scanning geometry. Therefore, one positive aspect of the Planck scan strategy is that one 
would not expect to see similar systematics coupling to the quadrupole for the Planck scan 
strategy. 

An an aside, our toy model can also be useful in evaluating scan parameters of future 
CMB satellites. For example, the Experimental Probe of Inflationary Cosmology {EPIC) [l2| 
mission proposal requires that pixel scans are as isotropic as possible to minimize polarization 
angle errors. We find that a spin and precession period similar to WMAP^ but with = 
45° and ^= 45°, would result in improved angular coverage at each pixel compared to 
WMAP Q. 

IV. OBSERVATIONS OF POINT SOURCES 

So far we have focused on the induced signal from the large CMB dipole component. 
Observations of point sources could also leave an interesting signature, depending on the 
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FIG. 3: Modulus of the scan tangent vector expectation value for Planck (left) and the induced 
map (right, units are mK) in Ecliptic co-ordinates. 




FIG. 4: Observed point source for the WMAP scan pattern, originally at an Ecliptic pole (left) 
and on the equator (right). We have set the angular pointing error to be larger than the beam size 
to exaggerate the effect. 

instrument beam size compared to the angular pointing error 6^. In Fig. [4] we show an 
example of two point sources, one originally located at an Ecliptic pole and one on the 
equator. We use the WMAP scan strategy and set the pointing error to be larger than the 
beam size to highlight the effect. At the pole, the source is observed as a circular ring with 
radius 6^^ while on the equator there are two arcs in the north/south directions. This is 
consistent with our earlier discussion of the direction of pixel scans at these locations. 

In reality, the beam size of WMAP is larger than the claimed pointing offset (for W-band 
the beam is ~ 12'), so one would observe a smearing of sources rather a distinct offset. We 
expect this effect would still be detectable since WMAP has a large catalog of ^ 400 sources 
at a signal-to-noise exceeding 5 Binning sources by Ecliptic latitude and then stacking 
could lead to tight constraints on timing residuals; however, we leave this study for future 
work. 

Point sources might also leave additional features in the maps if a timing residual was 
present. When operating in 'total intensity mode', we find no other features in the map apart 
from those localized around the point source, but in 'difference mode' we find low amplitude 
stripes in the map at a distance from the source corresponding to twice the beam opening 
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angle. These features are only present when there is a timing residual, and presumedly 
arise from the inability of the maximum likelihood solver to recover a single solution which 
is actually the difference between two maps. One could imagine in principle adding an 
additional free parameter to the map-maker to solve for the timing offset. However, since 
the gross effects of the claimed timing mismatch are not apparent in the data, there seems 
little point in investigating the more subtle effects. 

Finally, we should mention a related issue, which is that a timing offset would be at least 
partially calibrated out in practice. The WMAP beam boresight is calibrated on twice- 
yearly observations of Jupiter. Therefore, the timing offset would be degenerate with the 
offset of a given detector in the focal plane (relative to the boresight) in that particular scan 
direction. Timing offsets in other scan directions would then only be partially corrected 
by the Jupiter calibration. This could be easily be accounted for in our code, however due 
to uncertainty in the precise details of the calibration (e.g. Jupiter is observed for 45 day 
intervals twice- yearly) , we have not explored this correction in more depth. 

V. CONCLUSIONS 

We have investigated recent claims in the literature that the WMAP quadrupole is sys- 
tematic in origin, arising from an offset between the recording of pointing and temperature 
data from the satellite. Due to the WMAP scan pattern coupling with the direction of the 
dipole field, this effect results in a strong Y2-1 mode in Ecliptic co-ordinates, which happens 
to be similar in phase to the actual quadrupole. We find that the size of the effect for an 
offset of 25.6 ms cannot be large enough to match the observed quadrupole. We have also 
described reasons why the claimed effect is unlikely to be in the WMAP data, in addition 
to the WMAP team insisting that any error in their timing could not be nearly as large 
as claimed. Nevertheless, because of the similarity between the induced and primordial 
quadrupole signals, one should exercise increased caution in interpreting the amplitude of 
the primordial component. Planck has a significantly different scan strategy and will soon 
provide an independent measurement of the quadrupole. 
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